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I. INTRODUCTION 

There being so many excellent publications that had 
covered the study of gravitational lensing, our justifi- 
cation for another general article on the subject comes 
from the fact that we present new expressions for the 
optical scalars and deflection angle for a wide variety 
of matter distributions in terms of the matter compo- 
nents. 

Gravitational lensing has became a significant tool 
to make progress in our knowledge on the matter con- 
tent of our Universe. In particular, there is a large 
number of works that use gravitational lensing tech- 
niques in order to know how much mass are in galaxies 
or clusters of galaxies. One of the most exiting results 
was to reaffirm the need for some kind of dark matter, 
that appears to interact with the barionic matter only 
through gravitation. 

The question in which there is yet not general agree- 
ment is on the very nature of this dark matter. The 
most common conception is that it is based on coUi- 
sionless particles WeiOS], and where the pressures are 
negligible. However in the context of cosmological 
studies, one often recurs to models of dark matter in 
terms of scalar fields |MGN00l rULL02l IALS031 INSSOlj . 
There is also the possibility that dark matter were de- 
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scribed in terms of spinor fields [GuOTj . 

One method to study the nature of dark matter con- 
sists in observing the deformation of images of galax- 
ies behind a matter distribution that is the source of 
a gravitational lens. 

The fact that gravitational lensing can be useful 
for the study of the nature of dark matter has been 
emphasized many times, in particular in respect to 
the question of its equation of state |FV06[ IBK03[ 
INMCCMIO . 

In many astrophysical situations, the gravitational 
effects on light rays is weak, and the source and ob- 
server are far away from the lens, therefore they are 
studied under the formalism of weak field and thin 
gravitational lenses. The basic and familiar variables 
in this discussion are shown in figure [T] 




FIG. 1. This graph shows the basic and familiar angular 
variables in terms of a simple flat background geometry. 
The letter s denote sources, the letter I denotes lens and 
the observer is assumed to be situated at the apex of the 
rays. 



In this framework the lens equation reads 

a.9 



(1) 



The differential of this equation can be written as 
dp^^A^Se^, (2) 
where the matric A^j, is in turn expressed by 

1 - K - 71 -72 



^ b — 



-72 1 - K -f 7i 



(3) 



where the optical scalars k, 71 and 72, are known 
as convergence k and shear components {71,72}, and 
have the information of distortion of the image of the 
source due to the lens effects. 
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It is somehow striking that in most astronomical 
works on gravitational lensing, it is assumed that the 
lens scalars and deflection angle, can be obtained from 
a Newtonian-like potential function. These expres- 
sions although are easy to use, have some limitations: 

• They neglect more general distribution of 
energy-momentum tensor Tat, in particular they 
only take into account the timelike component 
of this tensor. In this way they severely restrict 
the possible candidates to dark matter that can 
be studied with these expressions. 

• They are not expressed in terms of gauge invari- 
ant quantities. 

• Since these expressions are written in terms of a 
potential function, it is not easily seen how dif- 
ferent components of Tat contribute in the gen- 
eration of these images. 

Moreover, most of them assume from the beginning 
that thin lens is a good approximation. 

In other cases in which the thin lens approximation 
is not used |BBV10| . the results are presented in a way 
in which gauge invariance is not obvious; however see 
IFKNOObj . 

In this paper we extend the work appearing in 
standard references on gravitational lensing [SEF92[ 
ISSE94I IWammbI IBarlO| and present new expressions 
that do not suffer from the limitations mentioned 
above. In particular we present gauge invariant ex- 
pressions for the optical scalars and deflection angle 
for some general class of matter distributions. In this 
first work on the subject, we study weak field gravi- 
tational lensing over a flat background. 

In section |ll] we present the general setting, start- 
ing from the geodesic deviation equation, where we 
fix some of our notation, and obtain gauge invariant 
expressions for the lens scalars. We also present a 
method of approximation for solving the lens equa- 
tions, that can be applied to any order. In section 
|III| and |IV[ we will restrict the study to axially sym- 
metric lens, and we present expressions for the lens 
scalars and deflection angle in the thin lens case. In 
section \V\ we concentrate in the spherically symmet- 
ric case, and after a general study of this geometry, we 
obtain expressions in terms of the energy-momentum 
distribution of these optical quantities. We end with a 
summary flnal section |VI| and a couple of appendices. 



II. INTEGRATED EXPANSION AND SHEAR 

A. General equations: The geodesic deviation 
equation 

Let us consider the general case of a null geodesic 
starting from the position ps (source) and ending at 



Po (observer). Let us characterize the tangent vector 
as £ = ; so that 

f^'V^r = 0; (4) 

that is, A is an affine parameter. 

We can now consider also a continuous set of nearby 
null geodesies. This congruence of null geodesies can 
be constructed in the following way. Let S' be a two 
dimensional spacelike surface (the source image) such 
that the null vector i is orthogonal to S. Next we can 
generalize £ to be a vector fleld in the vicinity of the 
initial geodesic in the following way: let the function 
u be defined so that it is constant along the congru- 
ence of null geodesies emanating orthogonally to S 
and reaching the observing point po. Then, without 
loss of generality we can assume that 

L = Vau; (5) 

which implies that the congruence has zero twist. 

We can complete to a set of null tetrad, so that m° 
and are tangent to S. At other points m° is chosen 
so that it is tangent to the surfaces u =constant and 
A =constant. Then a deviation vector at the source 
image can be expressed by 

c;'^ = <jTO° -ffm'^. (6) 

In order to propagate this deviation vector along the 
null congruence one requires, that its Lie derivate van- 
ishes along the congruence; that is 

= 0; (7) 

which is equivalent to 

£''Wb'^'' - <r^Vhr = 0. (8) 

From this it can be proved that 

^(4^') = 0. (9) 

The expansion and shear of the congruence are 
defined [Pir65] respectively by 

- ^v.r (10) 

and 

W\ = y^V(A)V-^''-02. (11) 

whose relation to the spin coefficient [GHP 73) quanti- 
ties p and a is given by 

0=-liP + P) (12) 

and 

|ap=aa; (13) 
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with certain abuse of the notation on the first appear- 
ance of \a\. 

Let us now calculate the covariant derivative of 
in the direction of £, 



while 



^f-Vb^r" = ^(cr)?fi" + qtVbrfi" + c.c. 



.^''Vfer = ^ml'Vit + c.c. 



where c.c. means complex conjugate. 
Let us note that 



and 



f-Vfer = (-/?' + /3)r - am" - pffi" 



i^Xl^rrf- = (e - e)m" - Rn" - tT; 



(14) 
(15) 

(16) 
(17) 



where we are using the GHP |GHP73j notation for the 
spin coefficients. In our case one has 



K = 0: 



(18) 



since t is geodesic. Notice that the Lie derivative of 
vector m in the direction of I is 

m]'' = (p + e - e) + am" + - /3 - f') t. 

Then, since the Lie transport of m in the direction of 
t should not have any I component; due to the fact 
that m are always tangent to the surfaces u ^constant 
and A ^constant; one obtains that 



Therefore, from eq.Q one has 

^(<^)m°+<;(e-e)m'' 
+ z. {am" + pfh") + c.c. = 0; 

which implies 

= £(?) + <r(e -e)+^p + ^a. 



(20) 



(21) 



(22) 



Using the GHP notation one can write the previous 
equation as 

= l>(<r) + crp + <rcr; 

where E> is the well behaved derivation of type {1, 1} 
in the direction of £. In order to have simple relations 
in terms of coordinate derivatives in the direction of 
A, the complex phase of m and fh can be chosen so 
that e = 0; so that finally one has 



dX 



-<;p - ?cr. 



(23) 



We see then that p determines the instantaneous ex- 
pansion, and a determines the instantaneous shear of 
the congruence. 



Let us recall from the GHP equations |GHP73] that 



e{p) = + cro- + $ 



00 J 



and 



P 



a p 



one has 

where Q is given by 

$00 *o 



*o ^-0 



with 



and 



Defining X by 



(24) 



£ia)^{p + p)a + ^o- (25) 
Defining the matrix P from 
p a 



(26) 



(27) 



(28) 



*00 = -;^i?a6«', (29) 



(30) 



(31) 



the equation for can be written as 

e{X) = -PX; (32) 

so that 

ii£{X)) = -QX] (33) 
which it only involves curvature quantities. 



B. Approximation method for solving the 
geodesic deviation equation 

Although the last equation can be integrated nu- 
merically without problems; it is sometimes conve- 
nient to have at hand some method for approximated 
solutions. So, next we present an approximation 
scheme that it can be applied to any order one wishes 
to obtain; although we will concentrate on the lin- 
ear approximation since in weak field lens studies it 
is consistent to consider linear effects of the curvature 
on geodesic deviations. 

Let us first transform to a first order differential 
equation. Defining V to be 



V = 



dX 
dA ' 



(34) 
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and 



and that 



one obtains 



^(X) 



X 



dX 



X 



V 



dX \ -QX 



(35) 



AX; (36) 



with 



A = 



I 
/ ■ 



(37) 



which gives 



Equation ( 36 ) can be reexpressed in integral form, 

(38) 



X(A)=Xo+/ .4(A')X(A')dA 



One can define the sequence 



X2(A) =Xo+ / ^(A')Xi(A')dA'; 



(40) 



and so on. 

Assuming that Q is in some sense small, one expects 
that this sequence will converge and therefore provide 
for the solution. 

Let us observe that 

X2(A) =Xi(A) 

r\ (41) 

+ / A(A') / A(A")rfA"rfA'Xo; 
Working out each term, one can see that 

X3(A)=^-/i/A>"'^^"'^^' 

V -/a>'^a' 



A{\') A{\") ^ A' A" 



-Q" 
-Q' 



(42) 



where we are using the notation Q' ~ Q(X'). Similarly 
one has 



A' A" A'" 



-Q' 
Q'Q'" 



(43) 



and 



A' A" A"' A" 



Q"Q"" 
Q'Q'" 



(44) 



Xi(A) = Xo+/ yl(A')XodA', (39) Q 



So one can see that only at the fourth product of ma- 
trices ^'s one has complete second order of matrices 



Returning to the sequence, the third element in first 
order is given by 



X3(A)=Xo + ^' J^dA'Xo 



+ 



A ,.A' 



Ao -^Ao 
A <.A' <.A" 



Ao -^Aq -'Ao 



-Q" 
-Q 

-Q' 







, ) dX"dX'Xo (45) 
dX'" dX" dX'Xo. 



^o)i-lU'jx"-Xo)Q"dX'dy\^ 

I-/,^,(A'-Ao)Q'dA' J °' 



(46) 



where one can check that the second row is just the doing an integration by parts as 
derivative of the first row. 

Let us note that in this equation one has not yet 
determined whether the position designated by A is to 

the future or the past of the position designated by Ao; r^' ^ 

so that one can use this approximated expression for / / (A"-Ao)Q(A")dA"dA' 

both cases, keeping the same direction for the vector ^° ^° ^ (47) 



i. If one changes the direction of the vector £, then 
one has to take into account that V changes to —V. 



(A- A')(A'- Ao)Q(A')dA'. 



Ao 



In particular, it is easy to see that (461 is invariant 
under interchange of A — > — A and V ^ — V. 

Note also that the double integral that appear in 
the first row and second column can by written by In the following we will make use of this equality. 



5 



C. The integrated shear and expansion 

Now in order to integrate the geodesic deviation 
equation, we must choose the correct initial condi- 
tions. In the case of hght rays belonging to the past 
null cone of the observer and intersecting S at the 
source, this initial conditions are X = and V ^ 0; 
since one can think the beam, starts backwards in 
time from the observer position, and so initially has 
vanishing departure, but with nonzero expansion and 
shear. 

Therefore in the linear approximation one has 

-\o)Q'd\']v{Xo); ^^^^ 



and 



V(A) =y- - ^o)Q' dX'j V(Ao). (49) 

In these integrations Aq indicates the position at the 
observer and from now on. As will indicate the position 
at the source. 

We observe from the first expression, that if the 
metric were flat {Q ~ 0), in order to get a deviation 



vector constructed from Xi , defined as X evaluated at 
Xs = Xq + ds, one must choose as initial condition 



V(Ao) 



1 



1 



(A. - Ao) 



X{X, = Ao+4) = -rXi. (50) 



Let us remark that we have just fixed the scale of the 
afHne parameter A to coincide with the measure of 
spacelike distances. 

But in the case of the presence of a gravitational 
lens, if an observer sees an image of "size" A'o, which 
means Xq = dgVo (since actually what is observed is 
Vo = V(Ao) ), then it should be produced by a source 
of size Xs = X{Xs), as described by equation (481 and 
depicted in figure 

In order to simplify the notation, we set from now 



on Aq = and As = d^, then eq.(48) reduces to 



Jo 



X'{ds~ X')Q' dX' ] Xo. (51) 



Note also that although a generic complex displace- 
ment should be ^ = l^je*"^, for our purposes it is 
enough to consider a complex displacement of unit 
modulus; namely <; = e^^. Then from equation (51) 
one would have 



A'(4-A')$oo(A')dA' 



where one can see that for the flat case one has 
^(A, (^) — e**^. From this equation, it is also observed 
that the expansion is only governed by the integra- 
tion of $00 1 and that the shear is only given by the 
integration of ^l^o. 



D. Expressions for the lens optical scalar in 
terms of Weyl and Ricci curvature from geodesic 
deviation equation 

In order to compare with the standard representa- 
tion of the lens scalar we note that the original de- 

I 



ds Jo 



A'(ds-A')*o(A')rfA' 



(52) 



r 



viation vector in the source will be given by eq.(51) 
i.e. 



if we make the following decomposition into real and 
imaginary part. 



To — To/?, + i'To/, 

^-0 = *oj? + «*o/; 



we obtain from cq. ( 53 ) that 



(54) 
(55) 
(56) 



/ f'^^ A' (4 - A') . , 



= 1 



, 1^00 ^or) 



dx] w - r ^'^V ^'^ %idX',,n- 



(57) 



Note also that in principle the integration must be ton in its path from the source to observer. How- 
made through the actual geodesic followed by a pho- 
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Flat lightcone 




Real lightcone 



FIG. 2. An object of typical dimension ds?" it appears to 
the observer to have a "size" ds<;S- 



ever the last expressions are valid only in the limit 
where the linear approximation is valid. If one con- 
siders a linear perturbation from flat spacetime, then 
the curvature components <i>oo and would be al- 
ready of linear order. Then, in the context of weak 
field gravitational lensing, it is consistent to consider 
a null geodesic in flat spacetime; since the actual null 
geodesic can be thought as a null geodesic in flat 
spacetime plus some corrections of higher orders. 

We choose then, a null geodesic coming from a 
source located at a distance ds from the observer, 
and select a Cartesian coordinate system where this 
geodesic will propagate along the y negative direction. 
As we mentioned previously, one can actually inte- 
grate the equations either along the physical direction 
or one can integrate to the past from a null geodesic 
that starts at the observer position. We make this 
second choice. 

We also need a null tetrad {1°", m", m°, n^}, adapted 
to this geodesic, 

r =(-1,0,1,0), 

m" =-^(0,1,0,1), 

1 (58) 
m'^=-^(0,-z,0,l), 

n'' =^(-1,0,-1,0). 

Now, in order to compare with the usual expressions 
for the lens scalars k, 71 and 72, let us recall that they 
are defined via the relation eq.([2]); but since it is a 
linear relation, one can relate the deviation vectors by 
the same matrix, namely 



where {c*, c*} are the spatial vector associated with 
{^s, ^0} respectively. In this expression, it is needed to 
determine the meaning of the indices of the two 
dimensional space of the images. In order to observe 
the natural Cartesian orientation we identify the first 
component of the two dimensional space with the z 
one of the complete system, and the second compo- 
nent of the two dimensional space with the x one. 
We need then, to know the components of the spatial 
vectors generated by <^o and similarly by in a 
Cartesian like coordinate system. In the case of , it 
is given by 



?^ =<rom° + Com'' 

= ^(^o(0,-i,0,l)+s^o(0,i,0,l)) 

=^(0,i(Co-';o),0, (<;o-hCo)) ^^^^ 



V2 



(0, w,0,<ro_R) ; 



and a similar expression is obtained for . 



Therefore, by replacing into eq.(59), we obtain 



■^si? = (1 - K - TO-^O-R - 72 W, (61) 

^s/ = -72'^o_R + (l-K + 7i)w; (62) 



which by comparing with eq.(57), implies that 



Jo 



71 



A'(4-A')$oodA', 



72 A'(4-A')*o/o'A'. 

Jo 



(63) 
(64) 
(65) 



Let us emphasize that these expressions for the weak 
field lens quantities are explicitly gauge invariant, 
since they are given in terms of the curvature compo- 
nents, which are gauge invariant. This is in contrast to 
the usual treatment of weak field gravitational lensing 
found in the literature, which use for example equa- 
tion (2.17) of reference |SEF92) as the source for the 
calculation of the lens scalars. 

Note that the last two equations can be written as 



7i + ^72 



1 



dj ''^'^ 



X')%dX'. (66) 



(59) 



As a final comment to this section it is important to 
remark that these expressions are valid for any weak 
field gravitational lens on a perturbed flat spacetime, 
without restriction on the size of the lens compared 
with the other distances. However, if we make use 
of the hypothesis of thin lens, these equations can be 
further simplified, as we will show below. 
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III. THE AXIALLY SYMMETRIC LENS 
(INCLUDING THE SPHERICALLY 
SYMMETRIC CASE) 

When one observes an astrophysical system very 
often one needs to extract information of the bulk of 
the matter distribution, which normally involves to 
make some model assumptions on the nature of the 
distribution. So, very often one considers spherically 
symmetric models or the less restrictive case of axially 
symmetric distribution. In this later case the axis co- 
incide with the line passing through the central region 
of the distribution and the observer. 

In this section then we consider the case of an axi- 
ally symmetric gravitational lens without introducing 
further assumptions on the extend of the lens. Later 
we will consider thin lenses. 

Using the same setting as in the last section, one is 
studying the motion of a photon which travels along 
the negative y direction, with impact parameter J and 
angle § from the z axis. Then one notes that the 
component $00 is a spin zero real quantity, and it 
only depends on the (J, y) coordinates. While the 
component ^0 is a spin two complex quantity and it 
has the functional dependence 



^-0 = |*o|e 



(67) 



where the phase is gauge dependent. For reasons that 
will become more clear during the study of spheri- 
cally symmetric systems, we define the real quantities 
'^Po{J,y) from 



From this we deduce that the optical scalars have 
the following dependence 

A«( J) = ^ [ ^ A'(4 - A')$oo( J, A') dX', (69) 
ds Jo 



and 

1 

71 +«72 = 



A'(4-A')^o(J,A')dA'. (70) 



This invite us to also define the real quantity j{J,y) 
from 



7i + *72 
so that one simply has 



-76 



„2ii9. 



(71) 



1 



A'(4-A>o(J,A')dA'. (72) 



IV. THE THIN LENS APPROXIMATION 

A. The general case 

Now, we will consider the case of a lens whose size is 
small compared with the distances to the source and 



the observer. Let there be again a Cartesian coordi- 
nate system such that the lens can be thought to be 
localized around the plane y = 0- 

Then as it was indicated in the last section, in the 
linear approximation we can replace the actual null 
geodesic by one in a flat spacetime. Then, considering 
a null geodesic as in the previous section, coming from 
a source located at a distance d^ from the observer, 
and at a distance dis from the lens, coming parallel to 
the y axis, but in the negative direction; we will use 
J to represent the impact parameter and i9 to denote 
the angle of the trajectory as measured from the z 
axis, in the {z,x) plane. We choose the scale of the 
affine parameter A such that the geodesic is described 

by 



(x(A), 2/(A), z(A)) = (xq, a - d;, zo); 



(73) 



i.e., A = indicates the position of the observer, and 
X ^ ds the position of the source. 

Then if we represent generically by C each one of 
the scalars {$oOj^o}j that appear in the expressions 
for the lens scalars, we have that by doing an integra- 
tion by parts one obtains the relation 



X'{ds - X')C{X')dX' 



X'ids - X')CiX') 



{ds - 2X')C{X')dX' 



{ds - 2X')C{X')dX'; 



where 



C{X') 



C{X")dX". 



(74) 



(75) 



Then if we assume a thin lens, the scalars C, will 
be sharply peaked around X — di, where it is located, 
and C(A) can be approximated by 



C(A') 





c 



yx<di-s 

VA >di+S; 



(76) 



where 5 di, 6 <ti dis and 5 <C ds- Therefore, we 
obtain 



{ds - 2A')C'(A')dA' C {ds- 2X')dX' 

Jdi 

= Cdi{di-ds) 

^ -Cdidis] 



(77) 



where we have neglected terms of order 0{^). Ob- 
viously all this also holds for the particular case of a 
delta Dirac distribution for the curvature components; 
however our relaxed notion of thin lens is entirely ex- 
pressed by the ( 76 ) behavior. 
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Finally we conclude that in the thin lens approxima- 
tion, the expressions for the lens scalars are reduced 
to 



where 



didis » 



didis » 
71 + *72 = — ^*o, 



Jo 

*o = / ^„dX, 
Jo 



(78) 
(79) 



(80) 



are the projected curvature scalars along the line of 
sight. 

We again emphasize that these expressions for the 
lens scalars are explicitly gauge invariants. 



B. The axially symmetric case (which includes 
the spherically symmetric case) 

1. The lens scalars m terms of projected Ricci and Weyl 
Scalars 

For axially symmetric lens (and in fact spherically 
symmetric lens), the projected curvature scalars are 
given by 

$oo(J) = / ' $oo(A')dA', (81) 



where one can see that 



MJ) 



-2ii) 



*o(A')dA'. 



(82) 



(83) 



The reason for the minus sign choice is that in 
many common astrophysical situations one would find 

MJ) > 0. 



By replacing in eqs. (78) and (791, we obtain for 
the lens scalars 



K = -^'^oo{J). 



7i 



72 = 



d, 
disdi 

ds 
disdi 



MJ) cos(2,9), 
MJ) sin(2i?). 



which implies that 



disdi -J . 

7 = --^MJ)- 



(84) 
(85) 
(86) 

(87) 



This equations can be compared to those of reference 
[FKNOOgi] : where they use different notation but sim- 
ilar content. 



2. Deflection angle in terms of projected Ricci and Weyl 
Scalars 

We wish now to express the deflection angle in terms 
of the curvature scalars. 
From eq.([2]) we know that 



4 



6' 



dis da"- 



disdi da"- _ 
dxi ' 



where in the last equality we have used that in the 
thin lens approximation 



We define the components of a* = (a^, a^) as 

ia^) = a{J)Cj,^); (89) 

since as we mentioned above we are respecting the 
Cartesian orientation in the two dimensional space of 
the images. We then obtain that the shears and con- 
vergence can be written as. 



K = - 



71 

72 

Noting that 



1 disdi / da^ da^ 

2 ds \ dzQ dxQ 
1 disdi ( dfx^ do? 
1 ds \ dzo dxo 
disdi da^ disdi da^ 

ds dxo ds dzQ 



xo =Jsin(i?), 
zq =Jcos(i?); 



we obtain 



71 



1 disdi / da 
\dJ 



2~d 



1 disdi 



72 = 



2 ds 
1 disdi 
2^ 



cos(2i9) 
sin(2?9) 



iJ) 

J 

da 
dJ 
da 
dJ 



aiJ) 

J 
a{J) 



It is interesting to note that 
K — 7i cos(2i?) — 72 sin(2i9) 



didis a{J) ^ 



(90) 
(91) 
(92) 

(93) 

(94) 
(95) 
(96) 

(97) 



ds J 

from which, using eqs. ( 84p6 ), it is deduced that 

a(J) = J($oo(J)+^o(J)). (98) 

It is worthwhile to remark that this constitutes an 
equation for the bending angle expressed in terms of 
the gauge invariant curvature components in a very 
simple compact form. We do not have knowledge of a 
previous presentation of this equation. 

It is also important to emphasize that we have de- 
rived the expression for the deflection angle from the 
information contained in the calculation of the optical 
scalars, coming from the geodesic deviation equation. 
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Note that, at first sight, it seems that if we recon- 
struct the lens scalars using eqs. ([94]|96|, from expres- 
sion (98) for the deflection angle, one would obtain 
some condition on the bending angle when compared 
with eqs. ( 84|[86 l; however this is only an apparent 
inconsistence. Let us see this in more detail. To be- 
gin with, by replacing eq. (98) into eqs. ( 94p6 ), one 
obtains 



disdi r 



2d., 



2{<Poo{J) + MJ) 



+J- 



dJ 



(99) 



71 



72 



disdi 
~2d7 



cos(2i9) J 




disdi . d^oo 
„ , sm 2p J — — 
2d, ^ ' \ dJ 




(100) 



(101) 



Proceeding with the calculation, we now use one 
of the Bianchi identities, as expressed in the GHP 
formalism [GHP 73] . namely 



l)*i - g'*o + g$oo - l>$oi = 0; 



(102) 



where in our case E> = l°'da, 3 = m'^da, and 3' 
fh°'da, i.e. 



d 








. d 




dx 


dz) 


. d 


d 


dx 


dz 



(103) 



that is, we are using the flat null tetrad system. The 
expression for the edth operator is correct due to the 
fact that the intrinsic two dimensional metric in the 
space (a;, z) is constant and therefore 3 = m. 

Let us now change to a polar coordinate system in 
the two dimensional subspace, so that 



d 


dJ d 


dd d 


dx 


dx dJ 


dx dS ' 


d 


dJ d 


dd d 


dz 


~lhdJ ~ 


dz dd' 



(104) 



with 



dd _cos{d) 

dx J 

dd _ sm{d) 



(105) 



In this case the metric of the two dimensional space 
(J, 1?) is not constant, so that in principle the edth 



operator acting on a quantity / of type (p, q) should 
be[GHP73] 

3/ = rn{f) + {~p(3 + qP')f- (106) 

but a direct calculation in the {J,d) frame gives all 
spin coefflcients zero. Therefore in this frame we also 
have 3/ = m{f ). Then we get 



1 

■7=2' 
1 

'-T2 



,ii5 



d 

dJ^ 
d_ 
dJ 



i d 
Jdd 
i d 
' Jdd 



(107) 



If we now project the Bianchi identity on the line 
of sight direction, i.e. by integrating along the y- 
direction, we obtain 



*ilo= -3'*o + 3$oo-*oilo^ 
and $01 ~ 



which, assuming "^i 
the lens, it implies 

3'(Vioe2^^) 
From this one finds 

,1pQ 



5(*oo). 



= 0; (108) 
far away from 

(109) 



dJ 



2^ = 



00 



(110) 



J dJ ■ 

Then by replacing this relation into eqs. ( 99|l01 ) we 
obtain eqs. ( [84p6 ), as anticipated. 

The Bianchi identities have not been used very often 
in the context of gravitational lenses; however we note 
that in reference [KK05, KC08 they have used them 
to obtain a Poisson like equation in order to determine 
the matter distribution. 

For the study of the errors committed in the use of 
the thin lens approximation one can read |KF08] . 



V. DETAILED STUDY OF STATIONARY 
SPHERICALLY SYMMETRIC LENSES 

Up to now we have presented gauge invariant ex- 
pressions for the deviation angle and the optical 
scalars in terms of the curvature components of the 
null tetrad adapted to the motion of the photons. 
In order to obtain expressions that use information 
of the structure of the sources one has to work with 
frames adapted to the geometry of the matter distri- 
bution which forms the gravitational lens. Therefore 
in this section we study the case of stationary spheri- 
cally symmetric sources. 



A. Spacetime geometry in standard coordinate 
system 

1. The metric 

For stationary spherically symmetric spacetime, the 
line element can be expressed by 

ds^ = a{r)dt^ - b{r)dr^ - r^{de^ + sin^ Odip'^). (Ill) 
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It is convenient to define <I>(r) and M{r) from 



a{r) 



„2<[.(r) 



and 



6(r) 



1 - 



2M(r) 



(112) 



(113) 



The more general distribution of energy-momentum 
compatible with spherical symmetry is described by 
an energy-momentum tensor given by 



T . ^ • 

(^1 2M(.)y 



Tee = Pt 



T^^ = Pt sin(0)^ 



(114) 



(115) 



(116) 



(117) 



where we have introduced the notion of radial compo- 
nent Pr and tangential component Pt- 
The Einstein field equations 



Gab = —8nTab, 
in terms of the previous variables are 



dM 
dr 



Attt g, 



(118) 



(119) 



where r* is chosen so that u is null. Then by inspec- 



tion of equation ( 111 ) one can see that 

dr* 



-dr; 



(124) 



du — dt —dr = dt — 

dr 

since then one has 
ds'^ = adu^ + 2^/abdudr-r'^{de'^ + sva^ edif^). (125) 

It is natural to define the principal null direction (.p 
from 



Ip = du; 

which implies that the vector is 



tp = g^'dub = 
where we have used that 

^gab^ ^ 2 d d Id d 



ah du dr h dr dr 

_]_,d_d_ _l d_d_ 

~r^^d6d6 ^ sm^ ed^d^ 

Let us define the null tetrad 

d 



: A 



dr' 



(126) 



(127) 



(128) 



(129) 



, d$ M + Airr^Pr 



dr 



1 - 



2M(r) 



d$ \ 2M 

^2d^ M dM 

dr r dr 

dM „ 
r-—+M^ Snr^Pt. 
dr 



The conservation equation is 

dPr , ,d<^ 2, 



(120) 



(121) 



(122) 



B. Geometry with respect to a null system 

1. The tetrad 

For our purpose, it is more convenient to use a null 
coordinate system to describe the spherically symmet- 
ric geometry. Let us introduce then, a function 



u = t — r* 



np = — -t- UA—, 
ou or 

with the complex null vector 

V2Po d 
r dC 



mp 



in terms of the stereographic coordinate C,. 
Therefore, one has 



1 



ah 



and 



U 



1 



26A2 2 ■ 



(130) 



(131) 



(132) 



(133) 



It is worthwhile to note that we have chosen to keep 
using r as a coordinate; which measures the surfaces of 
the symmetric spheres. Instead one could have chosen 
to use an affine coordinate f so that one would have 
£ = ; but then the surfaces of the symmetric spheres 



(123) would be some function of f , different from Airr^ . 
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C. The spin coefficients scalars and curvature 
components 

For the spherically symmetric metric, the non van- 
ishing spin coefScients are 



(134) 
(135) 
(136) 

(137) 
(138) 



p = 


1 

r 








UA 
r 










dPo 


^ .dPo 


V2r ( 


dy2 




fi' = 


^/2r ( 


dPo 
dv2 


dys 


i' ^ 


1 , dV 

-A — : 








2 dr ' 







where we are using for the stereographic coordinate, 
the decomposition ( = ^{y2 + Ws)- 

The curvature components that are different from 
zero are: 



$11 



$00 = -- 



A dA 
dr ' 



1 dAdU , 1„, ,,,, . , 



1 dU 

— i 

2r dr 



$22 = 



AC/2 dA 
r dr ' 



(139) 



(140) 



(141) 



Also let us note that 
du 

which then implies that 
1 d 



d_ 
Ft' 



Tip = 



2dt 



fa d_ 
b ^dr' 



(146) 



(147) 



In these last equations (^)|^ is meant at constant 

t; as opposite to the previous equations in which 
was meant at constant u. 



D. Spinor Ricci components in terms of 
energy-momentum components in the 
non-isotropic case 

The spinor Ricci components can be written in 
terms of the energy-momentum distribution as 



(148) 

(149) 
(150) 
(151) 



These expressions are exact for the spherically sym- 
metric spacetime. If one needs linear expressions 
around flat spacetime, one must set a = 1. 
Note that one has 



$00 = 


47r , 

— U? 

a 


+ Pr), 


$11 = 


TT{g- 


Pr + 2Pt) , 


$22 = 


an {g - 


^Pr), 


A = 




Pr~2Pt); 



$2 



C/2$i 



00- 



(152) 



and 



1 dAdU^ 
12 drdr^ 
1 /IdA 
3~d^ 



^VHU)A' 



AU - 



1 dU ^ 
6 dr 



12' 



^'2 = - 



1 dAdU 
6 dr dr 
1 /IdA 



A - -V2(C/)A2 
6 



2dU 



(142) 



(143) 



Note that from (124) and (1261 one has that 



and therefore 



ip = dt- 



a at 



-dr. 



ah ^ dr^ 



(144) 



(145) 



Using the expressions for $ii and A one can prove 
that 

- AdiUA) 1 , , 1, 

+ +2) (153) 

^2TT{g~Pr). 

Also, from the relation of the null tetrad compo- 
nents with the old variables, one can obtain that 



UA^ 



M{r) 



(154) 



This equation gives U in terms of A and M. 
Using this in the expression for $00 one obtains 



1 dA _ Anrig + Pr) 
A ~ 



(155) 



This is a useful equation only involving A; which al- 
lows it's calculation in terms of the components of the 
energy-momentum tensor. 

The contracted Bianchi identity (2.37) of |GHP73| 
for spherically symmetric metrics is 

I>$ii+D''loo+3I>A= (p'+^')$oo+2(p+y5)'lii. (156) 
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or explicitely, 



— ^+AC/— ^ + 2— A$oo + 3 — 
ar ar ar dr 

= -2 $oo-4-$ii, 

r r 



(157) 



which gives the conservation equation in the form 
ar 

where we are using 



(159) 



This is a very simple relation for ^2('') in terms of the 
energy density £'(r), the spacehke components and the 
mass function M{r). Our expression generalizes those 
of reference ' Dye77| for the case of anisotropic energy- 
momentum tensor. 



F. The bending angle and lens scalars in terms 
of energy-momentum components, curvature 
components and M(r) 



1. Relation between the scalars curvatures in the two 
different tetrads 



E. Simple relation for Weyl component »I'2 

Let us observe that 

*2 + 2A= 160 

r dr 



Then from equation (1531 one can deduce that 



47r, , M 



(161) 



In order to express the function a{J) in terms of the 
curvature scalars defined with the spherically symmet- 
ric tetrad, we need to know how the tetrads transform 
between them. To do so, let us recall that at lin- 
ear order, we only need the transformation between 
the flat tetrad {/°, to"^, to", rt°} adapted to the null 
geodesic coming from the source, and a flat tetrad 
{l'^ , m"^ , rh"' , n°'} obtained from {lp,Thp,fhp,np} by 
setting a = b = 1. Then, using standard spherical 
coordinates we have 



r = (1,-,^,-) ^ (l,sin(6')cos((/)),sin(6')sin(0), 008(6*)), 
r r r 

n° = (1, - sin(6') cos(0), - sin(6') sin((/)), - cos(6')), 

m° = Bail"") = ^(0, - cos(0) cos{(l>) + ism{(l>), - cos(6') sin(</)) - z cos((/)), sin(^)), 
V 2 



(162) 
(163) 
(164) 



^Boin = ^(0,-cos(6l)cos(<?!))-isin(<?!)),-cos(6')sin((/))-f icos((/)),sin(6')). (165) 
V 2 



In these expressions we use the symbols 3o and 3o to 
denote the edths operators of the spheres of symmetry, 
with unit radius. 

The transformation tetrad will be of the form 



Cimw" -I- c„-n". 



— Qn^" ~ C;^TO 
TO — Cjjinl ^mrh^ Cjnrh^ ^mf^ ' 



TO 



(166) 
(167) 
(168) 
(169) 



where the notation is cin — l°'na, and so on. From 
these relations we can construct the transformation of 
the Ricci and Weyl scalars, but it is more convenient 
and easy to work with the spinor diad associated to 



r 



the tetrad. Then the transformation of the dyads will 
be 



(170) 



together to the condition that o^, conform a spino- 
rial base, i.e., 

AD-BC^l, (171) 
From the relations given in Appendix A, we get 



A =4=v/l - sin(6l) sin(,^)e*'^ 
V 2 

B =v/l + sin(6i) sin((?!))e' 



C =-v/l + sin(6') sin((/))e*' 



(172) 



D = — ^l- sin(6l) sin(0)e" 
V 2 
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where rj and rj' satisfies 



-cos(6l)sin(</)) +icos((/)) 
e = . (173) 

^l~sm^{9) sin2(,/>) 
^ cos(6>) +^sin(6l) cos((/)) ^^^^^ 
^/l -sin2(e') sm^{(j)) 
z + ix 



J 



= cos(i?) + i sin(i^), 



and in the last equaUty it was used the fact that J = 
r-y/ 1 — sin^(6') sm^{(j)), (see Appendix A). Note then 
that, V = -d. 

The general transformation between tetrads induces 
the following transformation on the curvature scalar 
$00 and *o, 

=A'^A^^oa + 2^2^5*01 + ^^5^*02 

+ 2AA^B^iQ + iAABMii + 2ABB'^^i2 
+ A'^B^^-,0 + 2B^BA^2i + B'^B'^^22, 

(175) 

=A^^o + 4:A^^i+6A'^B^i'2 (176) 

In the spherically symmetric case these transforma- 
tions simplify considerably and finally, at linear order 
one has $22 = j^qo so that 



vfo = 3^^2(r)e2'^ 



2,P ~ 1 ~ 
^00 = ^(*ii - 7*00) + 



00- 



(177) 



(178) 



Note that in this case, the integration is on the coor- 
dinate y, instead of using arbitrary afBne parameter. 
Also note that r — \/ + y^- 

This constitutes an important explicit relation for 
the bending angle in terms of the curvature as seen in 
an spherically symmetric frame; which is the natural 
frame for the sources of the gravitational lens. 

3. Expressions for the bending angle in terms of 
energy-momentum components and M[r) 



Using eqs. pl8| , pl9| , ( [I6I] ) and p79| ) we get 
an expression for the bending angle in terms of the 
mass, energy density, and spacelike components of the 
energy-momentum tensor, namely 



a{J) = J 



-d, 



3J2 fM{r) 



in 



(180) 



+4TT{g{r) + Prir))]dy. 



This is a new and useful relation for the deflec- 
tion angle in terms of the physical fields which are 
the sources of the gravitational lens. It is also worth 
mentioning that this expression for the bending angle 
can also be deduced from the geodesic equation using 
standard techniques, as it is shown in Appendix B. 

It is curious that the bending angle does not depend 
explicitly on the tangential spacelike components of 
the energy-momentum tensor. 



4- The optical scalars in terms of spherically symmetric 
components of the curvature 



From equation (|84|) and ( 87 1 one obtains 
2J2 



and 



disdi 



-($1 



7^00) 



disdi 



3^^2(0 



dy. 



dy, 
(181) 

(182) 



2. The deflection angle in terms of spherically 
symmetric components of the curvature 



These expressions give the optical scalars in terms 
of gauge invariant expressions for the curvature com- 
ponents adapted to the symmetry of the matter dis- 
tribution which is the source of the gravitational lens. 



From (98), the function a(J) expressed in terms of 



the spherically symmetric null tetrad reads, 



a( J) = J 



-7^*2 + ^($11 -4 



*oo dy 



^^00) 



(179) 



5. Expressions for the lens scalars in terms of 
energy-momentum components and M (r) 

In a similar way, the lens scalars, in terms of the 
spherically symmetric physical fields, are given by 
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K = 







ds J 


-d, 


didu 




ds J-d 





P + Pr 
3M 



Pt - Pr) 



dy 



Anip + Pt- Pr) 



(183) 



dy. 



These new expressions let us see explicitly the con- 
tributions of different components of the energy- 
momentum tensor on the optical scalars. One can 
see that a couple of terms disappear in the isotropic 
case in which Pr = Pt- 

Our expressions are valid for generic energy- 
momentum distributions, but it is usually not con- 
sidered in the literature the possible implications of 
non vanishing spacelike components of Tab- In fu- 
ture works, we will consider the implications of mod- 
els with non-trivial energy-momentum tensors on ob- 
served gravitational lenses. 



G. Two simple examples 

In order to show the application of our treatment of 
gravitational lens, we will consider next two standard 
models that are often used in representing the source 
of gravitational lenses. 



1- A monopole mass (Schwarzschild) 

As a simple example let there be a monopole distri- 
bution characterized by a mass M, therefore a simple 
computation gives I'oo = 0, and ^2 = — then by 
considering that the observer and the source are far 
away, one can replace in the extremes of the integra- 
tion (as is usually made) — >■ 00 and d/ — > cxi, then 



$00 

■00 



0, 



3J2 ~ , 4M 

-^^2dy=^, 



(184) 
(185) 



ily obtain the well known results 

AM 

K = 0, 

didis 4M 



and by replacing into eqs.(84), (87) and (98), we read 
1 kno 

a{J) 



7 



J2 



(186) 
(187) 
(188) 



2- The isothermal profile 

One simple model of dark matter that is used to ex- 
plain the rotation curves of galaxies is the isothermal 
profile, which is defined by the density function 



P = 



(189) 



where Vc is the circular velocity. 

Since Vc <C c, the pressures in this model are negli- 
gible. Then we obtain. 



00 



jdy 



J 



-dy 



r** " J 

From these relations follow the well know results, 

a — 27rvl, 

didis vItt 

ds J 
didis Tivl 



(190) 
(191) 



(192) 
(193) 



7 = 



J 



(194) 



VI. FINAL COMMENTS 



Several works on gravitational lensing reach up to 
the expressions that relate the optical scalars with the 
curvature components in terms of the tetrad adapted 
to the motion of the photons; we have here also pre- 
sented expressions for the bending angle in terms of 
the curvature components. Furthermore, we have pre- 
sented above expressions for the optical scalars and 
deflection angle directly in terms of the matter com- 
ponents of the sources of the gravitational lens, valid 
for an extended class of matter distributions. In order 
to do that one has to assume some structure for the 
source, so that in this first work on the subject, we 
have treated the first natural model of spherical sym- 
metry for the sources. But in section III we have pre- 



sented expressions that are valid also for spheroidal 
distributions, since we only required axissymmetry 
along the line of sight. 

Our expressions circumvent several deficiencies as 
are: gauge dependence, lack of explicit expres- 
sions, neglect of spacelike components of the energy- 
momentum tensor, etc. It is probably worthwhile to 
remark that since the function M{r) is determined in 
terms of the Q{r) by equation (119), all our expres- 



sions are explicit expressions in terms of the energy- 
momentum components of the matter generating the 
gravitational lens. As a trivial check of our equa- 
tions we have presented two simple examples for which 
the optical scalars and deflection angle are readily ob- 
tained. 

The extension of this study to sources with different 
structure and to the cosmological background will be 
presented elsewhere. 
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Appendix A 
1. Tetrad transformations 



First, we note that 



J2 J2 



J2 



r2 J2 j^y2 j2 _|_ J.2 sin(0)2 sin((?!))2 ' 
and solving this for J'^ /r"^ we find 
J2 



(Al) 



= 1 - sin(6')2 sin(0)2 



(A2) 



The complete equations that satisfy the spinorial com- 
ponents are given by 

AA = cm = ^ (-1 + sin(0) sin(0)) , (A3) 



AB = -ci,^ ^ ^i'^^'^) - i cos((/.)) , (A4) 



BB = c„~ = (-1 - sin(6') sin(0)) , (A5) 



CC = c„ft = -(-l-sin(0)sin((/))), (A6) 



rip) — _„ - _ 



(cos(6') sin((/)) - i cos{(j))) , (A7) 



^C* — c„,,, — 



2^2 



(— cos(0) — i sin(6') cos((/))) , 



(A8) 



AD = -c^ 



BC = -c„ 



1 



1 



Appendix B 

1. Deflection angle in terms of Tab from geodesic 
equation 

The four velocity vector of the particle has modulus 



1 



1 



,dr , 



2M y J \ ' 



(Bl) 



where A is an affine parameter of the geodesic, and we 
have already made use of the symmetry that allows us 
to study just the motion in the equatorial plane 9 — ^. 
The constant k has values 1 for massive particles and 
for massless particles. This choice for k sets the 
unit for the affine parameter for the massive particle 
case; however the unit for the massless case remains 
undetermined. 

There are also two integrals of motion. J is a con- 
stant of motion associated to the existence of a rota- 
tional Killing vector which can be expressed by 



J 



(B2) 



E is another constant of motion associated to the ex- 
istence of a timelikc Killing vector, which can be ex- 
pressed by 



E 



2<S> 



dt 
dX' 



Then equation (Bl ) takes the form 



or 



(d^^2 

^^'^dx' 



J2 



(B3) 



(B4) 



e-'^E^ 



(1- 



2M, 



k(1 



2M. 



(B5) 



(sin(6') -I- sin((/)) — i cos{9) cos((/))) , which can also be expressed as: 



(A9) 



(sin(6') — sin(0) — icos(6') cos(0)) , 

(AlO) 



dX r 



2M 



£;2e-2*(i 



2M^ 
r 



(B6) 



BD = c, 



ml 



'V2 



{cos{e) + ism{0)cos{(t>)), (All) 



DD = c^i = ^ (-1 + MO) sin(0)) ; (A12) 
and its complex conjugates together to the condition 
AD-BC=1. (A13) 



It is observed that the choice of the affine parameter 
A is related to the definitions of the constants of mo- 
tion J and E. Since tends to zero in the asymptotic 
region, it is natural to take A so that E ~ 1. This is 
equivalent to say that in the asymptotic region one 
has dt = dX. 

In this way there is no more freedom in the choice 
of units for J. For an incident photon traveling in the 
—y direction, with coordinate x = xq, the Newtonian 
expression for the angular momentum, for a unit mass 
particle gives J = rvsm{ip + f ) = rcos{(p) ~ xq; 
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that is with this choice of afBne parameter, J has the 
meaning of asymptotic impact parameter xq. 

For convenience in the algebraic manipulation, let 
us define 



ai($) = l- 



-2$. 



(B7) 



so that e 2* = 1 — ai{^) in the above equation. 
Therefore, for a photon, one can express (B6) by 



,dr,, J2 „2M 2M ( 2M\ , , 

+ ^ - + +1 "1 * ^ 1- 

d\ r'' r \ r ) 

(B8) 

The corresponding potential for the motion of a 
photon is 



Vi. = - J- 



MM,, , , 
ai($); (B9) 



y 



r dr 



(B14) 



with the velocity notation — ^ &n dvy = ^a- 

Let us assume the initial conditions: x = xq^ y » 
2M , fx = and Vy = —1. Then, in this case, xq is 
the impact parameter, so that J = xq- 

After passing through the gravitational lens, the 
trajectory will be deflected so that, in the asymp- 
totic region one would have Vx\oo = Sv and Vy\oo — 
— — (Sv)^; since the photon must travel at the ve- 
locity of light. Then the bending angle can be calcu- 
lated from 



a — — arctan ■ 



-Sv. 



(B15) 



which, we remark, is an exact expression. 

If one considers only linear departures from the flat 
metric, one whould replace e~^* sa (1 — 2<i>); and so 
one would obtain 



+2 



kM 
r 



M 



= E^ - K. 

(BIO) 



For the massless case, by choosing the parametriza- 
' 1 for r — )■ 00 one has E = 1, and 



tion A so that || 
therefore, one defines 



J^M M 



(Bll) 



The motion of a photon can be deduced from the 
Lagrangian 



L 



dr , 
dX' 
.dx , 
d\' 



(B12) 



where in the last equality we have used Cartesian like 
coordinate system with r = \pxP^^\^]p- . This system 
obviously has the integral of motion ( B2 1 with La- 
ir 

grangian energy c — 



The equations of motion are: 



dvx 



X dVf 
r dr 



(B13) 



The variation in the velocity can be calculated from 



dr 



:dVi 

dr 



dy; (B16) 



where we have taken dX = dy. 

The coordinate system has origin at the center of 
the spherical symmetry. In the approximation of a 
lens contained in a plane, the center is in this plane. 

To consider the equation of motion of a massless 
particle in the more general case, we also use the equa- 
tions of motion in the Cartesian like coordinate sys- 
tem, where now the potential is given by (Bll); so 
that 



dr 



J2 M 



J^M M 



J2 



dM 
dr 



J2 



M 



dr 



r^ 
M 



+ {l ^)47rrp 



(B17) 



iirr-^Pr 



and in the linear regime one has 



dVi 
dr 



r 



M (r) 47r 



Anr {p{r) 



3 

Prir)) 



p{r) 



(B18) 



Finally the deflection angle is given by 



a 



Sv 



-d. 



xo dVe , 
r dr 



Xq 



3J2 
r 



M (r) 47r 



+ -7rP{r) + 47rr (p(r) -|- Pr{r)) 



(B19) 



dy] 
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where it is understood that r = 

This coincide with expression ( 180 ) appearing 
above. 
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